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Abstract. We describe a particular aspect of the effects of the parent cluster tidal 
field (CTF) on stellar orbits inside cluster Elliptical galaxies (Es). In particular we dis- 
cuss, with the aid of a simple numerical model, the possibility that collisionless stellar 
evaporation from elliptical galaxies is an effective mechanism for the production of the 
recently discovered intracluster stellar populations (ISP). A preliminary investigation, 
based on very idealized galaxy density profiles (Ferrers density distributions), showed 
that over an Hubble time, the amount of stars lost by a representative galaxy may sum 
up to the 10% of the initial galaxy mass, a fraction in interesting agreement with ob- 
servational data. The effectiveness of this mechanism is due to the fact that the galaxy 
oscillation periods near equilibrium configurations in the CTF are comparable to stel- 
lar orbital times in the external galaxy regions. Here we extend our previous study to 
more realistic galaxy density profiles, in particular by adopting a triaxial Hernquist 
model. 



1 Introduction 

Observational evidences of an Intracluster Stellar Population (ISP) are mainly 
based on the identification of intergalactic planetary nebulae and red giant stars 
(see, e.g., [1],[2],[3],[4],[5]). Overall, the data suggest that approximately 10% 
(or even more) of the stellar mass of clusters is contributed by the ISP [6]. 
The usual scenario assumed to explain the finding above is that gravitational 
interactions between cluster galaxies, and interactions between the galaxies with 
the gravitational field of the cluster, lead to a substantial stripping of stars from 
the galaxies themselves. 

Here, supported by a curious coincidence, namely by the fact that the char- 
acteristic times of oscillation of a galaxy around its equilibrium position in the 
cluster tidal field ( CTF) are of the same order of magnitude of the stellar orbital 
periods in the external part of the galaxy itself, we explore the effects of interac- 
tion between stellar orbits inside the galaxies and the CTF. In fact, based on the 
observational evidence that the major axis of cluster Es seems to be preferen- 
tially oriented toward the cluster center, N-body simulations showed that model 
galaxies tend to align, as observed, reacting to the CTF as rigid bodies [7] . 
By assuming this idealized scenario, a stability analysis then showed that this 
configuration is of stable equilibrium, and allowed to calculate the oscillation 
periods in the linearized regime [8]. In particular, oscillations around two sta- 
ble equilibrium configurations have been considered, namely: 1) when the center 
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of mass of the galaxy is at rest at center of a triaxial cluster, and the galaxy 
inertia ellipsoid is aligned with the CTF principal directions, and 2) when the 
galaxy center of mass is placed on a circular orbit in a spherical cluster, and the 
galaxy major axis points toward the galaxy center while the galaxy minor axis 
is perpendicular to the orbital plane. 

Here, prompted by these observational and theoretical considerations, we 
extend a very preliminary study of the problem [9], by evolving stellar orbits 
in a more realistic galaxy density profile: for simplicity we restrict to case 1) 
above, while the full exploration of the parameter space, together with a complete 
discussion of case 2), will be given elsewhere [10]. It is clear, however, that both 
cases are rather exceptional. Most cluster galaxies neither rest in the cluster 
center nor move on circular orbits, but they move on elongated orbits with 
very different pericentric and apocentric distances from the cluster's center; in 
a triaxial cluster many orbits are boxes and some orbits can be chaotic. These 
latter cases can be properly investigated only by direct numerical simulation of 
the stellar motions inside the galaxies, coupled with the numerical integration 
of the equations of the motion of the galaxies themselves. 



2 The physical background 

Without loss of generality we assume that in the (inertial) Cartesian coordinate 
system C, with the origin on the cluster center, the CTF tensor T is in diagonal 
form, with components Tj (i = 1,2,3). By using three successive, counterclock- 
wise rotations (ip around x axis, $ around y' axis and tp around z" axis), the 
linearized equations of the motion for the galaxy near the equilibrium configu- 
ration can be written as 

.. AT 32 AI 32 .. AT 31 AI 31 •• AT 21 AI 21 

ip= ip, •&= i9, ip= V, (1) 

H H ^3 

where AT is the antisymmetric tensor of components ATij = Ti — Tj, and Jj are 
the principal components of the galaxy inertia tensor. In addition, let us also 
assume that T x > T 2 > T 3 and h < h < h, i.e., that AT 32 ,AT 31 and AT 21 
are all less or equal to zero (see, e.g., [8], [10]). Thus, the equilibrium position 
associated with cq. (1) is linearly stable, and its solution is 

<p = (fM cos(ui,pt), d — $m cos(w#t), ip = ipM cos(atyt), (2) 

where 



AT 23 AI 32 /AT 13 AI 31 /AT 12 AI 21 



— h — ' ^ = V — T 2 — ' ^ = V — T 3 — ' (3) 

For computational reasons the best reference system in which calculate stellar 
orbits is the (non inertial) reference system C in which the galaxy is at rest, 
and its inertia tensor is in diagonal form. The equation of the motion for a star 
in C is 

x' =K T x-2f2 Av' - SI Ax'- f2 A {fi Ax 1 ), (4) 
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where x = 1Z((p, "4>)x', and 

fl = (ip cos ■d cos tp + flsinip, — ^> cost? sin f/> + i? cos V 1 , ¥> sm- $ + V 1 )- (5) 
In eq. (4) 

"K T ai = -Va,'^g + (n T TTZ)x', (6) 

where 4> g (x') is the galactic gravitational potential, W x > is the gradient operator 
in C", and we used the tidal approximation to obtain the star acceleration due 
to the cluster gravitational field. 

3 Galaxy and cluster models 

For simplicity we assume that the galaxy and cluster densities are stratified on 
homcoids. In particular, the galaxy density belongs to ellipsoidal generalization 
of the widely used 7-models ([11], [12]): 

< \ - M s 3 ~ 1 1 m 



where M g is the total mass of the galaxy, < 7 < 3 and 

= 22^-, ai>a 2 >a 3 . (8) 
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The inertia tensor components of a generic homeoidal density distribution 
(in the natural reference system adopted in eq. (8)), are given by 

h = — aia 2 a 3 (a? + al)h s , (9) 

where h g = p g (m)m 4 dm, and so I\ < I2 < Note that, from eqs. (3) and 
(9) it results that the frequencies for homeoidal stratifications do not depend on 
the specific density distribution assumed, but only on the quantities (ai, cti-, olz)- 
We also introduce the two ellipticities 

— = l-e, — = l-n, (10) 
ai ai 

where e < rj < 0.7. 

A rough estimate of characteristic stellar orbital times inside m is given by 
Porb( m ) — ^Pdyn( m ) = V '37T / 1 Gp~^(m) , where Pg( m ) is the mean galaxy density 
inside m. We thus obtain 



P orb (m) ~ 9.35 x 10 6 J — ^ m^ 2 (l + m)( 3 ^/ 2 yrs, (11) 
V M sM 

where M g; n is the galaxy mass normalized to 10 n M Q , is the galaxy "core" 
major axis in kpc units (for the spherically symmetric 7 = 1 Hcrnquist model [13], 
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Fig. 1. Distribution of the dmax/di ratio vs. di/cei after an Hubble time for the model 
galaxy at rest, di is the initial distance of the star from the galaxy center, while d max 
is the maximum distance from the galaxy center reached during the simulation. 



R e ~ 1.8c*i); thus, in the outskirts of normal galaxies orbital times well exceed 
10 s or even 10 9 yrs. For the cluster density profile we assume 

^ c ( m ) = 7T~r^2V2 ' ( 12 ) 
(1 + m^Y 

where m is given by an identity similar to eq. (8), with a\ > ai > 03, and, in 
analogy with eq. (10) we define 02/01 = 1 — fx and 03/01 = 1 — with fi < v < 1. 
It can be shown (see, e.g., [8], [10]) that the CTF components at the center of a 
non-singular homeoidal distribution are given by 

Ti = -2TiGp C fiW l {^, u), (13) 

where the dimcnsionless quantities Wi are independent of the specific density 
profile, wi < W2 < W3 for a\ > 0,2 > 03, and so the conditions for stable equi- 
librium in eq. (1) are fulfilled ([8], [10]). The quantity p Cj o is not a well measured 
quantity in real clusters, and for its determination we use the virial theorem, 
M c <Ty = — U, where <j\ is the virial velocity dispersion, that we assume to be 
estimated by the observed velocity dispersion of galaxies in the cluster. Thus, 
we can now compare the galactic oscillation periods: 



2tt _ 8.58 x 10 8 o 1|25 q 

Uip yj(y — pi){j] — e) cry, 1000 



2ir 8.58 x 10 s a 1|250 
Pff = — ^ = ■ yrs , (14b) 

U)$ y/PT] (Ty, 1000 

2vr 8.58 x 10 8 ffll|250 
Pi, = — ^ == ! yrs . (14c) 
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(for small galaxy and cluster flattenings, where 01,250 = ai/250 kpc and cry. 1000 = 
cry/lO 3 km/s, [10]) with the characteristic orbital times in galaxies. Thus, from 
eq. (11) and eqs. (14abc), it follows that in the outer halo of giant Es, stellar 
orbital times can be of the same order of magnitude as the oscillatory periods of 
the galaxies themselves near their equilibrium position in the CTF. For example, 
in a relatively small galaxy of M gi n = 0.1 and a± t x = 1, P or b — 1 Gyr at m ~ 10 
(i.e., at ~ 5i? c ), while for a galaxy with M gi n = 1 and a±,i = 3 the same orbital 
time characterizes m ~ 7 (i.e., ~ 3.5i? ). 




Fig. 2. Distribution of the d max /rfi ratio vs. dk/ai after an Hubble time for the same 
galaxy model as in Fig. 1, when oscillating around its equilibrium position in the CTF. 

In order to understand the effects of the galaxy oscillations on the stellar 
orbits, we performed a set of Monte-Carlo simulations, in which we followed 
the evolution of 10 4 - 10 5 "1-body problems" over the Hubble time by inte- 
grating numerically eq. (4). At variance with [9], where we used simple and 
easy-to-integrate Ferrers density profiles, here we study orbital evolution in a 
more realistic (but also more demanding from the numerical point of view) 
galaxy density profile, namely a triaxial Hernquist model, obtained by assuming 
7 = 1 in eq. (7). The gravitational potential inside the galaxy density distribu- 
tion, in a form suitable for the numerical integration, was obtained by using an 
expansion technique useful in case of small density flattenings ([10], [14]). The 
initial conditions are generated by using the Von Neumann rejection method 
in phase-space (for details see [10]): note that, at variance with the analy- 
sis [9], now "stars" are characterized by initial velocities that can be different 
from zero. The code, a double-precision fortran code based on a Runge-Kutta 
scheme, runs on GRAVITOR, the Geneva Observatory 132 processors Beowulf 
cluster (http: / /obswww.unige.ch/~pfennige/gravitor/gravitor _e.html). 
The computation of 10 orbits usually requires 2 hours when using 10 nodes. 
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4 Preliminary results and conclusions 

We show here, as an illustrative case, the behavior of the ratio d max /d{ as a 
function of d;/ai, for a moderately flattened galaxy model (e ~ 0.2 and r\ ~ 0.3), 
with M g = 10 11 M Q , semi-mayor axis cvi = 3 kpc, and maximum oscillation 
angles equals to 0.1 rad. The cluster parameters are ai,250 = cv,iooo = 1, A* = 0.2, 
v = 0.4, and the total number of explored orbits is N to t — 10 4 . In order to show 
the effect of oscillations, in the following simulations we artificially eliminated 
the direct contribution of the CTF, as given by the second term in the r.h.s. of 

' <■:• (6). _ 

In Fig. 1 we show the result of a first simulation in which the galaxy is not 
oscillating: obviously, the ratio d max /di is in general (slightly) larger than unity, 
due to the initial velocity of each star. In Fig. 2 we show the result for the same 
galaxy model, when oscillating around the equilibrium position: the effects of 
the galaxy oscillations are clearly visible as a global "expansion" of the galaxy. 
As a reference, the solid line indicates the expansion ratio required to reach the 
representative distance of 10i? o from the galaxy center. Thus, it is clear that the 
galaxy oscillations are certainly able to substantially modify the galaxy density 
profile. In particular, it will be of interest the study of the (more realistic) case in 
which the galaxy is in rotation around the cluster center. In this case we expect 
a different behavior of stellar orbits as a function of the distance of the galaxy 
center of mass from the cluster center: in fact, while inside the cluster core the 
CTF is compressive (see, e.g., [7], [8]), outside the CTF is expansive along the 
cluster radial direction, and in this latter case its direct effect should increase 
the expansive effect due to the galaxy oscillations. These cases are discussed in 
detail in [10]. 
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